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Calibration of microwave noise sources requires a standard source, an attenuator, and 
an instrument for comparing microwave power levels. A modified Dicke radiometer is 
used for this comparison. f~A detailed analysis of the radiometer is undertaken employing 
the cascading matrix of each waveguide element, in order to evaluate the balance error. 
The analysis accounts for the general case of noise power arising from lossy elements of the 
radiometer and shows how these can often be accounted for by discussing temperatures in 
excess of some reference temperature T . It also displays several potential sources of error 
in such radiometers and provides the quantitative basis for keeping these sufficiently small. 
Fluctuations are analyzed as a basis for estimating sensitivity of measurement. Procedures 
for adjusting the radiometer to account for the characteristics of the two input paths are 
given. The results are carried through second order in the various reflection coefficients or 
related scattering matrix elements involved, which are assumed to be small quantities. ~"\ 

1. Introduction 

The National Bureau of Standards is concerned with the development and maintenance of 
national standards of microwave noise power and with calibration of working or interlaboratory 
standards. This work requires a standard source, a source to be calibrated, an attenuator, and 
a radiometer for detecting the equality of microwave noise power levels [Dicke, 1940; Selove, 
1954; Mayer, McCullough, and Sloanaker, 1958; Estin, Trcmbath, Wells, and Daywitt, 1960; 
Wells, Daywitt, and Miller, 1962]. One source is attenuated a known amount until its power 
level is equal to the power level of the second source. \\\ order to evaluate the calibration, all 
errors introduced by the radiometer must be known to second order in the small error quantities. 
No such analysis is found in the literature. 

The radiometer is first described. An analysis of the waveguide portion is undertaken 
using the cascading and scattering matrices of network components. Noise arising in the net- 
works is explicitly included. The mixer-IF is analyzed, and spectra at various points in the 
receiver are described. The effects of statistical fluctuations are estimated. The analysis 
yields a balance equation which is specialized to three convenient methods of operation. A 
typical error analysis is performed with particular scattering matrices. Finally, all errors are 
evaluated numerically and summarized. 

2. Description of Radiometer 

One source of microwave noise power is an argon gas discharge tube located within a 
rectangular waveguide [Mumford, 1949; Johnson and Deremer, 1951; Knol, 1951; Bridges, 
1954; Knol, 1957]. The noise power produced by a gas discharge may be calculated if thermo- 
dynamic equilibrium is assumed; however, such an assumption is not justified. Such calculated 
results could be in error by about 0.5 db. The noise power from such a device is stable, however, 
to better than 0.001 db when the discharge current is regulated. The argon noise source 
employed for WR 90 waveguide (8.20 Gc/s to 12.4 Gc/s) consists of a tube about 1 cm in diameter 
and 30 cm long mounted at a small angle to the axis of a rectangular waveguide, intersecting the 
broad wall of the waveguide. 

The most suitable standard of microwave noise power is a blackbody radiator [Sutcliffe, 
1 956], the power of which is given by the low frequency approximation to the Planck radiation 
law. One configuration may consist of a uniform piece of sintered silicon carbide mounted in a 
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waveguide section, as shown in figure 1, and heated to a uniform temperature of about 1,000 °K. 
Water cooling is used on one end of a short stainless steel thermal insulating section, so the 
temperature drop is taken across a short distance (1 cm). 



PLUG 
WAVEGUIDE 
- LOSSY TERMINATION 



Figure 1. A typical standard may consist of a lossy 

termination in low-loss waveguide. 

The standard is heated in an oven. 
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The radiometer is patterned after that of Dicke, with ferrite modulators in place of his 
rotating vane attenuators (figs. 2a and 2b). In critical applications, motor-driven rotary vane 
attenuators may serve as modulators, as ferrite modulators may have drifting characteristics 
and are sensitive to stray magnetic fields. Each input is alternately connected to the hybrid 
junction, after which the signal from either source travels the same path. The modulated 
noise signals are heterodyned to an intermediate frequency by a local oscillator and amplified. 
The signal is further detected, filtered by a narrow-band filter centered at the modulation fre- 
quency, amplified, and ultimately displayed on a recording galvanometer. 
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Figure 2a. Generator portion of a radiometer. 
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Figure 2b. Receiver portion of radiometer. 



(b) 

We shall derive an adjustment procedure that compensates many of the inequalities in the 
two branches of the radiometer. The analysis is done in detail so as to clearly display all the 
assumptions and approximations that are used. Other waveguide configurations are easily 
analyzed by the matrix methods presented here. 

A calibrated attenuator is used to attenuate the source with higher power until the unknown 
source and the known source balance each other. The final calibration is in terms of decibels 
above JcT Av [IRE, 1960], where k is Boltzmann's constant, T is an accepted reference temper- 
ature, and Av is the limiting RF bandwidth of the System. The maximum resolution of source 
temperature is limited b}^ system noise and the bandwidths of the system, which determine 
output fluctuations. 

3. Radiometer Analysis — Waveguide Portion 

The radiometer waveguide network shown in figure 2a is a member of a general class of 
waveguide networks which consists of the connection to a common junction of arbitrary groups 
of series-connected two-port components. Each group may be represented by its equivalent 
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cascading matrix, which is the product of the cascading matrices characterizing each two-port 
component. Such a network is shown in figure 3a. We shall be interested in the voltage at 
reference plane A-A. We shall therefore represent each cascading matrix and source by an 
equivalent source of reflection coefficient Tt, and equivalent temperature T' i} as shown in figure 
3b. For consideration of the voltage at reference plane A-A, the network to the left of the 
reference plane may be replaced by an equivalent source of impedance z G , reflection coefficient 
T G , and containing an equivalent generator of voltage e G , as shown in figure 3c. 
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Figure 3a. General n-port network with losses, 
equivalent to figure 2a. 

Each of n arms contains a total equivalent cascading matrix P t at a 
temperature T\ The junction is at a temperature Tj. We shall 
be concerned with the voltage at reference plane A-A, shown in 
arm n. 



Figure 3b. This general n-port network is the one 
in which the sources and cascading matrices shown 
in figure 8a have been replaced by equivalent termi- 
nations at temperatures T/, source voltages e/, and 
with reflection coefficients Y/ . 

The reference plane A-A is shown. 
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Figure 3c. In this figure, the entire network of 
figure 3b to the left of reference plane A-A has been 
replaced by an equivalent generator of voltage e G > 
with an equivalent source reflection coefficient Tq 
which is related to the equivalent generator impedance 
z G . 

The waves a' and &« are discussed in the text. The direction of these 
waves is intentionally chosen as shown in order to correspond to 
the incident, and emergent waves at the reference planes A-A in 
figures 3a and 3h. 



Ill figure 2b the mixer-IF output voltage is postulated to be a linearly related to the voltage 
at its input terminals. The rest of the receiver may operate further on this voltage 4 to give an 
output indication related to this voltage, such as the modulation amplitude or its mean square, 
apart from superimposed statistical fluctuations, of course. 

We must calculate the total terminal voltage at reference plane A-A as shown in figure 3. 
We shall solve for this voltage in figure 3c, then generalize to figure 3b, and finally generalize 
to figure 3a. This solution may then be used with appropriate simplifications to analyze the 
particular configuration of figure 2a. 

At any kih reference plane of the network, the real, time dependent fields, if restricted to 
a finite time interval often conveniently normalized to unity, may be represented by a Fourier 
integral. Provided the interval is sufficiently long, the Fourier amplitudes become known in mag- 
nitude in the sense that they have a spectral density (on the axis of frequency v) relatable to kT. 
Each Fourier component of the fields of a given mode may, from well-known waveguide theory, 
be discussed by means of amplitude coefficients v k , ik, Uk, and b k . These coefficients are closely 
related to voltages, currents, and incident and emergent voltage waves, respectively. These 
quantities, applicable to a single assumed propagating mode, are related through the charac- 
teristic impedance z (here assumed the same for each arm) by the equations 



v=a+b, 
z i=a—b, 



(3.1) 
(3.2) 



where the subscript k has been dropped. 



141 



We may therefore discuss transmission through the network using complex amplitudes of 
harmonically time dependent voltage density waves at frequency v and synthesize to the actual 
case of nonsinusoidal time dependence of noise by the usual integration over the appropriate 
frequency range. The network transmission properties will be characterized by matrices whose 
elements are complex functions of v. The mean square voltages over the time interval or the 
powers may be obtained by ParsevaFs theorem, that is, an integration of the square of the 
absolute value of the (normalized) Fourier amplitude over a frequency range which spans the 
limiting bandwidth of the system. This relation makes it convenient to deal with rms rather 
than peak values of a, b, and v in what follows. 

The waves a 1 and b l in figure 3c may be decomposed into two parts, given by 

a t =a+a°=T L b t +a°, (3.3) 

b^b+b^T^+b . (3.4) 

The quantity a or b° is the wave that would be launched down a reflectionless line by the partic- 
ular generator, L or G, respectively. For a , a reflectionless line would mean T G —(z G — z Q )/ 
(z<?+2o)=0, or z G =z , and similarly for b°. Then 

a°=e L ZoK2 Q +z L )=i(l^T L )e L} (3.5) 

b =e Q Zo/(z +z G )=%(l-T G )e G . (3.6) 

It may easily be shown, using (3.1), (3.2), and (3.4) that 

v=a i +b t =(l + T G )a Q /(l-T L T G ) + (l + T L )bV(l-T L T G ). (3.7) 

T G and T L are independent of frequency for sufficiently small Av. T L is fixed for a particular 
receiver input impedance, and T G will not vary during the process of making a measurement by 
adjusting the calibrated attenuator if the radiometer is properly designed. These observations 
permit us to absorb the coefficients of a and 6° in (3.7) into the constants of the receiver re- 
sponse and to focus attention on a and b° only. To the extent that the equivalent reflection 
coefficient 1^ of the source is not constant as the radiometer is operated, an error is incurred. 
If the variation of T G is known, a correction could be applied. If a is constant also, the receiver 
output is, apart from a constant, proportional to b°. 

This situation is generalized to the general n-port network of figure 3b. The equivalent 
source e G results from that portion of the network to the left of the reference plane A-A in 
figure 3b. Let there exist at each port of the network an incident wave a 1 and an emergent 
wave b l . We index the ports by i=l, 2, . . ., n. 

Consider the total incident wave a\ at port i. By the superposition theorem it consists 
of the sum of contributions from each source of the network operating alone. We index the 
sources by X=0, 1,2,..., n, where the index zero designates sources in the junction, and the 
other indices designate the sources terminating the respective arms. For \^i, the contribu- 
tion a,i\ is that partial wave incident on the iih reference plane which has been established by 
the source X^i. This partial wave, incidentally, is a reflection from T-. For \=i, it is con- 
venient to further decompose the partial wave established by the source X=i into parts a\ and 
a H where a t is the wave that would be launched by the source i down a reflectionless line, and 
a u is the remainder of the actual partial wave established by the source X=i. Only the portion 
a u is a reflection from r •. These considerations may be expressed as 

x=o 

The total emergent wave b\ at port i likewise consists of the superposition of contributions 
from all sources in the network. For X^O, the contribution b iX is that wave emergent from 
the iih reference plane which has been established by the source Xf^O. This partial wave, 
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incidentally, is a scattering of all the incident waves. For A=0 it is convenient to decompose 
the contribution from the sources within the junction into parts b\ and b i0 , where b\ is the wave 
which would be launched by the junction source (A = 0) down arm i if T- = for all i, and b i0 
is the remainder of the actual partial wave established in arm i by the source 0. Only the 
portion b m is a scattering of incident waves. Thus the expression 

b\=V\+ib &<A=ft?+»i (3.9) 

represents the total emergent wave at port i. 

The analysis preserves the familiar reflection law for the partial waves in the form 

a iX =T f i b i x ) X^O, 

a i o=T , i (b° i + bio), X=0. 

Summing over all A we have, using (3.8) and (3.9), 

at^T'iW+bi). (3.10) 

The analysis also preserves the familiar scattering law for the partial waves, for scattering 
matrix elements s ik , in the form 

bix=J2 s ifc a k x+s iX (a x +ax^), A^O 

n 
k = l 

bi0~2_j S ik a k0> A = 0. 

Again summing over all A we have 

6<=2 s ik (a° k +a k ). (3.11) 

/< = 1 

In matrix notation (3.8), (3.9), (3.10), and (3.11) become 

a l =a°+a (3.12) 

b l = W+b (3.13) 

a=T'b l (3.14) 

b=Sa t , (3.15) 

where r' is the diagonal matrix with elements T f t , i=l, 2, . . . n; S is the nXn scattering 
matrix; and the a's and 6's are column matrices. These are the basic equations of the circuit 
and constitute the generalization of (3.3) and (3.4). 

These four equations are easily solved for the matrix b t , eliminating the unknown matrix a. 
We obtain, using the unit matrix 1, 

6 f -(l-Sr / )- 1 (Sa°+& ). (3.16) 

From the matrix form of (3.1) applied to a 1 and b l , the voltage at any port is given by 

v=a +(l + T')b t =a +(l + T')(l-ST')- 1 (Sa° + b ). (3.17) 

This is the general solution for an w-port network. The voltage density at any port 
contains contributions from the terminations of all ports and from the losses within the junction 
itself. This shows the possibility of "leakage" from terminations other than the desired 
one if these terminations are not sufficiently isolated. 
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Now we may further generalize the situation to that of figure 3a, where scattering matrices 
have been introduced into each arm which will be used to perform the modulation and the 
balancing. We have defined the scattering matrix by (3.15). For a two-port network we 
may define a cascading matrix P such that 



\bj \p 2 i P22) \a 2 ) 
It is related to the scattering matrix by 

\P2l V2lJ \«11 



I/S21 — S 22 /S2l\ 

/*a s 12 — S11822/S21/ 



(3.18) 



(3.19) 



This matrix has the property that the cascading matrix for a number of two-ports in series is 
the product of their respective cascading matrices, allowing us to deal with any combination 
of microwave components in each arm, even including lossy sections of line. The total cas- 
cading matrix may then be changed back into the total scattering matrix by the relation 

All *12\ /2>2l/Pll P22—pV2P2l/pu\ , g ^ 

\*2i Sn) \ 1/pn —P12/P1J 

We wish to know the equivalent source reflection coefficient V and the equivalent source 
voltage density e' of a two-port driven by a source of reflection coefficient r and of voltage 
density e. (The index of the two-port has been temporarily suppressed.) It can easily be 
shown that 

r / =s 22 +s l2 s 21 v/(i-s 11 r). (3.21) 



One can also show that, neglecting sources within the two-port matrix, 

e'=s 21 (i-r)e/(i-6< n r)(i-r'). 

The matrix a has elements of the form of (3.5). Thus 

a?=i (i-r{)«{=J 4(i-i\W(i-shr,). 



(3.22) 



(3.23) 



The elements s l 21 are explicitly displayed since the comparison of power levels described later 
will depend on these elements. The matrix 6° will depend on the properties of the junction 
and any source in the junction. So long as these properties are time-independent, the matrix &<> 
need not be explicitly evaluated for our purposes. 

Having developed this very general formalism, we now apply it to the three-port junction 
shown in figure 4. In this case (3.16) becomes 



U= T 31 a?+ 2*3208+ Tnal+Nnbl+NvM+NzM 



(3.24) 
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Figure 4. General three-port network consisting of 
three equivalent scattering matrices S 1 , S 2 , and S^, 
of temperatures T 1 , T 2 , and T 3 connected to a com- 
mon junction characterized by a scattering matrix S 
(without superscript) and of temperature Tj. 
We wish to calculate the voltage at reference plane A-A . 
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where the coefficients are given by 

T zl =N zl s n +N Z2 s 21 +N zz s zl (3.25) 

7 32 = N zl s 12 + N Z2 s 22 + N sz s, 2 (3 .26) 

T zz =N zl s lz +N Z2 s 2Z +N zz s Z2 (3.27) 

iv 81 =[«»«nrlri+« a r 1 / (i-* la ri)]/det (i-ST') (3.28) 

iNT n =[*u««iriri+*tari(i-* 11 rO]/det (i-Sr r ) (3.29) 

iv 33 = [-*u« M rir;+ (1 -%r0 (1 -^rOJ/det (1 -Sr') (3.30) 
det (i-j8T0 = (i-«nr0(i-*ara(i— *^^^ 

-sia^rlr^i-^rO-^^rar^i-snrO-s^irir^i-sssrs). (3.31) 

We have previously noted that 1^ is constant. If also, through careful design and use of 
isolators, T[ and T f 2 are independent of changes in the quantities of (3.21) caused by modu- 
lation or substitution of different terminations r, all 7 and N coefficients in (3.25) through 

(3.31) are constant; and, if the //,' are constant, then 

H= 7 31 a?+ 7 32 ^+ 7 33 aS+constant. 

If further a% is constant, as it would be if the elements of S z are constant and contain no vary- 
ing sources, then 

H= 7 31 a?+ 7 32 a£+ constant. (3.32) 

For the radiometer problem, all sources in the network consist of thermal noise, so that 
we shall always be interested in the absolute squares of the waves, which are relatable to tem- 
peratures. Furthermore, waves arising from different thermal sources, such as a? and a° in 

(3.32) will be uncorrelated so that the absolute square averaged over any finite band will con- 
tain no contribution from cross products. With the understanding that this averaging will 
always eventually be done, such cross products will not be carried hereafter. 

The spectral density of voltage for grey-body noise is [Allis and Herlin, 1952] 

|<f=2jfc7Re (2)=2m (i-|r! 2 )/|i-ri 2 , (3.33) 

for real z } We use the factor 2kT instead of the more commonly seen factor ikT because 
we extend the density over the axis of negative as well as positive frequency. Thus in figure 
3c the absolute square launched waves are 

\a«\z=Ul-\T L \>)kT L z , 

m^ui-iTomToz,. 

Then in figure 3b the absolute square elements of the matrix a have the form 

KI 2 =Ki-!r;! 2 )&7;.2 . 



1 If the characteristic impedance is complex, then zo in (3.33) and following work must be replaced by 

zo->(Re zo)[l+2 Im z Im r/Re z (l-|r| 2 )] 

as may be verified by direct calculation of Re z in terms of r. The correction factor is very close to unity. If the correction is carried forward to 
(3.38) and (3.39), and thence to (6.1), it will almost completely cancel, except to the extent that differences in z and T i for arms 1 and 2 alter the 
slight departure from unity of the above expression. This correction is in practice smaller than the corrections retained and is therefore dropped 
al this point. 
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If the two-ports of figure 3a are either lossless or at °K, they will make no thermal 
contribution to e[ or T\ of figure 3b. The temperature T[ equivalent to e\ is obtained in 
terms of T t and the properties of the two-port from (3.22) and (3.33). The absolute square 
incident waves at the reference planes of the junction are then 

k?l 2 =il^il 2 (i-|r,| 2 )^,c /|i-^r,l 2 . (3.34) 

However, when the scattering matrices are each at some temperature T\ (3.22) must contain 
an additional term proportional to T\ Equation (3.22) so modified and (3.33) give 

I g «l- |il1i!r^i-r;p 2kT t z ^^ 2 +A l 2kT%=2kT' t z (i-\T^W-T' i \ 2 (3.35) 

where T t is the temperature of the ith source of reflection coefficient I\ T l is the tempera- 
ture of the 'ith scattering matrix, A t is a coefficient to be determined, and T[ is the temperature 
of the ?"th equivalent source of reflection coefficient 1^, as in figure 3b. We can evaluate the 
coefficient A t by noting that, when T l =T i =T Q . we must also have by the Second Law of 
Thermodynamics Ti=T . Dividing by 2kT Q z and rearranging, we obtain 

A= _i-|r;| 2 |«i|*(i-|r f |») 



|i-r;| 2 |i- Sl \r,.| 2 |i-r;| 2 

Substituting this value into (3.35) we may write 

T >jr, Kl 2 d-|r«| 2 ) -| ,, r Kl 2 d-lr f l 2 ) 1 T 
4 L ii— «fmr*|*Ci— |ril a )J "Tji-aMW-TOJ ' 

or, equivalently, 

(r '-^ =(r '-^[ ii-1i!w-i?;i') } ■ (3 - 36) 

which reduces to the former value of T\ for T*=0. 

Hence in figure 3a the absolute square incident waves are in general 

|o?| a =i ifcr^a— |r{|»)+i fc(T*— r*)zo(i— ir^i^i^ivii— *{ir f | a , (3.37) 

which are to be substituted into the absolute square of (3.32). The terms ^ |T r 31 | 2 (l — | rj| 2 )^rT 1 ^ (> 
and \\T^{X — \Y'^)kT 2 z^ are constant if T[, rj, T\ and T 2 are constant. 

For convenience we will absorb certain constants of proportionality into new definitions 
characteristic of the equivalent sources on arms 1 and 2 of figure 4: 

(3.38) 
and 

(3.39) 
Also, let 

and 



which are constant to the degree s\ x and s 2 n are constant. We shall also wish to compare the 
temperatures of the two sources mainly in terms of the p n elements of the cascading matrix 
defined in (3.19). Hence we define 

Pll=l/«21 

and 

2n=l/*2i. 
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We can now write the absolute square of (3.32) as 



i^i 2 =ir.;i 2 A-.+m' 2 - 1 ^ 



' bnl 2 



iQnl 



(3.40) 



apart from an additive constant. 

The important result is that, as modulation of p n and q u is introduced, the time dependent 
part of I&3I 2 depends on the excess temperatures {T l —T l ) and (T 2 —T 2 ). 

4. Radiometer Analysis — Receiver Portion 

We now describe the transmission of 63=63 + 63, emergent from the tee, through (he 
mixer-IF network. Physically, we must account for the following: possible image suppres- 
sion, frequency conversion of both signal and image, filtering effect of IF bandwidth, and 
addition of noise power by the mixer and IF. From the mathematical description of the 
superheterod3me process, the mixer should be considered a three-terminal-pair device, in 
which inputs at signal and image frequencies must he considered, as well as the output at 
the intermediate frequency. 

We formulate the system of equations describing (he mixer and IF as follows. Let e 
be the incident and / the emergent waves, functions of frequency; let the matrix U character- 
ize the scattering properties; and let the arguments v \ v be indicated by the subscripts ±, 
that is, e+ = e(v-\-vo) = e+(i>), e 2 = e 2 (v), etc. The mixer-IF network shown in figure 5 may be 
described by 
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(4.1) 
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FlGURE 5. Mixer 



IF network, described by tin 
matrix U. 






Matrix elements u 3i and u z2 are similarly functions of frequency, and may be written as the 
product of an IF amplitude response i{v) and an HF response r± = r(v± v )=r±(v). Thus 
u n =i(v)r + , and u 32 =i(v)r_. Amplitude responses are indicated by lower-case letters, and 
power response by capital letters, related thus: x(v)x*(v)=X(v), where * denotes complex con- 
jugate. Thermal noise has been treated by considering 63 and a- 6 ; all other noise is accounted 
for by the quantity n. For conventional mixers, Ui 2 =u 13 =u 2 i=u 23 =(). 

In the absence of any sources beyond the IF terminals we further define T by 6 2 =r/ 2 . 
Thus the wave emerging from the mixer-IF is 



f 2 =i(v)[r + bU + r-b^+n}/(l-Tu 33 ). 
We obtain using the rms amplitudes in (3.40) 



/.= 



1 — r^ 3 3\ 2hi+ 



Pi 



2n+ 



ffn- 



.„) 



+ constant. 



(4.2) 



Using Parseval's theorem, let us average the power incident on the mixer-IF output ter- 
mination over times long compared with the J F period but short compared with the modulation 
period by writing 

P=^ 1 [~ fijidr. (4.3) 
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We note that we are concerned only with sums of absolute squares, as a + , a_, c + , and c_ 
are all uncorrected. 

We note also that corresponding terms taken in + and — by pairs are symmetric about 
v=0, as are I(v), |1 — r^ 33 | 2 , and \n\ 2 , so that the integrand as a whole is symmetric in v, as it 
must be, and may be replaced by twice the integral from to °o . 

We also suppose that I{v) is relatively so narrow in the neighborhood of v if that all other 
quantities in the integrand may be replaced by their mean values oyer the band about v if . 

We also choose to assume local symmetry about v = v in the absolute squares of a, c, 
T'n, T' i2 , pn, and q n , considered individually as functions of the arguments i>±v . Specifically, 
we assume at v=v if) 

\a(v if ±vo)\ 2 =\a(—v if ±v )\ 2 . 

Then b} T virtue of the symmetry of \a\ 2 about the origin of its argument, imposed by the reality 
condition, namely, 

|a(, + ,o)| 2 =|a(-,-, )| 2 , 
we have 

\a + \ 2 =\a_\ 2 = a 2 

and similarly for the other named quantities. 2 

We do not assume such local symmetry in the KF response R, in order to account for any 
RF frequency dependences, such as introduced by possible image suppression filters perhaps 
placed in arm 3 of the junction. Nevertheless, the sum R+ J r R_ is symmetric about v=0. 

Applying these considerations to (4.3), but temporarily leaving Z? + + i?_ under the integral, 
we obtain 



p= ! 



Zo\l-Tu 33 \ 2 ' 



n 



' 12 



- M 2 + 
Pii\ 



£j32 
tflll 



)T U(v)d v + ,. * l2 \n\ 2 T /«^+constant, (4.4) 

J J -co ^Ol 1 — L U M\ J -co 

where the quantities outside the integral are evaluated at v=v if and where 

U(v) = [R(v+v ) + R(v-v )]I(v). (4.5) 



We now introduce the equivalent IF bandwidth referred to the response at v if , 

_J 

= 2/( 



VTT^ j Kv)dv. (4.(5) 



Finally, extracting J? + +7?_, we have 

P= 2I(R + +B^ V » + OTV 2/A, |«| 2+ con 8 tant, (4.7) 

O |i— r^ 33 | 2 V IphI 2 |gn| 2 / z \l-Tu zz \ 2 ' ' 

where the quantities are evaluated at v if . 

This is the power emergent from the IF amplifier in terms of the waves a and c leaving the 
two noise sources which are to be compared, and contains the effects of the waveguide attenu- 
ators and modulators and all noise. 

The following illustrates the power spectrum at various points in the receiver to provide 
a better qualitative understanding of receiver operation. The receiver may observe an RF 
spectral region U'(v), qualitatively like the rectangular pair example in figure 6. This result 
is available by simple and obvious changes of variables in the function U(v) in (4.5). 3 



2 If such symmetry does not exist in the spectral characteristics of a and c, then they must be treated as R immediately below. Lack of sym- 
metry would cause errors in the radiometry unless circumvented by image suppression. 

3 In the right-hand member of (4.5), first make the change of variable i»'=v-f-y in the term R(v+v ) 1(f), to obtain R(v')I(u'—v ). Then make the 
change v"=v— j> in the term R(v—v^I{y) to obtain R(v")I(v"~ v ). Dropping the primes and adding, we obtain the new function of p, 

U>(v)=R( v )[I{ v + v ,)+I( v -v,)]. 
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The power represented by (4.7) is slowly modulated at an audiofrequency v m through the 
time dependence of pn and q n . The IF amplifier is usually followed by an envelope detector. 
Near the origin the spectral distribution of the output of an unfiltered second detector of the 
quadratic, quadratic envelope, or linear envelope type is qualitatively similar to the example 
Q(v) in figure 7. Quantitative results for tin 4 quadratic envelope and linear envelope detectors 
are given in Law son and Uhlenbeck [1950). 

An audiofrequency filter at v m of width Av selects the fundamental modulated component 
of the second detector output voltage, designated by subscript m, and rejects the d-c and 
second harmonic components. The spectral density W(v) after the audiofrequency filter is 
qualitatively similar to the example shown in figure 8, where the frequency scale has been 
appropriately expanded. 

To obtain a d-c output to drive a chart recorder, the modulation filter output may be fur- 
ther rectified by a third detector. Analysis may be carried through quantitatively for the case of 
no modulation signal, depending on the particular type of detector. For either a quadratic en- 
velope or linear envelope detector, the spectral distribution W'(v) is qualitatively similar to 
figure 9, with details available in Lawson and Uhlenbeck [1950] as to the dependence of the 
spike area and triangle height on W(v m ) and B G . If a signal at v m is present due to the modu- 
lation, the spike area will greatly increase; and there will also be a continuous rectangular 
spectrum arising from heats of signal and noise, also described in Lawson and Uhlenbeck 
[1950]. To further reduce output fluctuations, the extent of the spectrum may be cut off at an 
arbitrarily small frequency (for example less than B Q ) by controlling the bandwidth of the low 
pass filter. 

The recorder is driven by this spectrum, the area under the zero frequency component 
providing the mean deflection and the nonzero components contributing to the residual fluctua- 
tions. The receiver thus finally indicates that Fourier amplitude at or near v m (or its square 
if there are two square-law devices) of the time-varying power given by (4.7). Specifically, 
the receiver is very insensitive to the constant terms in that equation. 
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Figure 6. Spectral region observed by receiver of rec- 
tangular pair response (no image suppression) . 
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Figure 7. Spectral density near the origin of output 
of square law second detector, fed by rectangular 
pair. 

If tlie input noise power is modulated corresponding to a tempera- 
ture change AT, signal peaks, qualitatively shown by the dotted 
lines, appear at the first and second harmonic, corresponding to 
a signal amplitude proportional to AT. A large d-c component 
exists [Lawson and Uhlenbeck, 1950]. The d-c component and 
the second harmonic peaks are discriminated against by a modula- 
tion frequency filter. This avoids practical problems of the usual 
d-c detection "such as d-c gain instabilities and could facilitate the 
obtaining of narrow post-detection bandwidths by synchronous 
detection at v m , for example. 
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Figure 8. Spectral density after the modulation 
frequency filter, G(V), of rectangular response. 
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Figure 9. Qualitative spectral density after a third 
detector fed by WW of figure 8. 

The analysis is that of a C\V signal in a hand of noise Apq [Lawson 
and CJhlenbeck, L950]. A large d-c component exists. Results 
are qualitatively similar for either a square-law detector or a linear 
detector. The dependence of tin 1 spike area and the continuous 
spectra on signal and noise may be found in the reference, for both 
cases. 
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5. Radiometer Output Fluctuations 

Output fluctuations should be analyzed to estimate to what extent, if any, they may limit 
the attainable accuracy. The general problem is that of modulated noise passing through 
successive stages of detection, any stage of which may consist of quadratic, linear envelope, 
or some other kind of detector. To get solutions for the purpose of either design evaluation 
or performance evaluation, the problem must be specialized to a particular configuration. 
Examination of many possible designs is beyond the intent of this paper; but for completeness 
a particularly simple example will be given, which is satisfactory for our purpose, and which 
is a slight extension of the usual analyses in that it discusses the modulation filter and the 
third detection. More refined analyses may be made using the appropriate references [Rice, 
1945; Lawson and Uhlenbeck, 1950; Colvin, 1959]. 

At the balance condition, even though |aj| 2 and |a°| 2 are modulated, \U\ 2 appears unmodu- 
lated. We shall calculate the fluctuations in the absence of modulation of |6|| 2 , and then 
compare them with a signal due only to the modulation of \bs\ 2 , assuming that the fluctuations 
do not significantly change when modulation is introduced. We shall assume that the input 
to the second detector is a noise amplitude, j(t), having Gaussian probability density, char- 
acterized by some noise figure N wdiich contains both the equivalent source temperature T G 
which produces 6|, and noise arising within the mixer — IF. We shall assume, in figure 2b, a 
quadratic second detector, followed by a modulation filter and by a linear third detector. 
The input to the square-law second detector consists solely of a noise voltage f(t) of spectral 
density F(v) given by the integrand of (4.3). The autocorrelation function of this voltage, 

^r=lim i f f(t)J(t+r)dt, 

7'->oo J- ,7 

is related to the spectral density by 

i/v-2 f F(v) cos 2irvrdv; (5.1) 

Jo 

and the mean square value of the voltage is 

^ =2 J F(v)dv. (5.2) 

Jo 

Rice [1945, his eq (4.10-1)] has shown that for a quadratic device [output q(t) 
=aj 2 (t)} the autocorrelation function of the output is given in terms of the autocorrelation 
function of the input bv 

*(r)=a 2 W Q +2tf). (5.3) 

The power spectrum is, since ^(r)=^(— r), 

Q(y)=2 [ ¥( r ) cos 2Tvrdr. (5.4) 

Conversely, since Q(v) = Q(—v), we have 

¥(r)=2 f Q(v) cos Znvrdv. (5.5) 

Jo 

We note in passing that practical square-law detectors responding to band-limited 
noise about some intermediate frequency do not often have this law, as they do not respond 
to the second harmonic component implicit in the law. Rather, they may be more properly 
called quadratic envelope detectors, having output proportional to the square of the envelope 
of /(£), that is, to the power averaged over many IF cycles. Nevertheless, it may be verified 
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that for band-limited noise such as from an IF amplifier, the output spectra of the quadratic 
detector and the quadratic envelope detector are identical except for contributions in the neigh- 
borhood of the second harmonic, which are filtered anyway. We may thus consider the 
analysis of the quadratic detector quite applicable. 

Let G(v)\ denote the power response of the modulation filter. Then the spectral density 
W(v) shown in figure 8 at the output of the modulation filter is 

W(v) = 6(v)Q(v). (5.6) 

Let us proceed to assemble W{v). From (5.4), (5.3), and (5.6) we obtain, after integration, 

W{v)=d 2 G{v)[tlb(v)+±^ Vr cos 2tvt(It], (5.7) 

using the representation 

5(j>)=2 j cos27nWr, (5.8) 

which has the property 

A(v)8(v)dv=A(0). (5.9) 



£ 



The convolution theorem states that the Fourier transform of the product of two functions 
is the convolution of their transforms. Thus, using * to mean convolute rather than complex 
conjugate, 

2f i/vAt cos 2ttpt(It= f F(v-p')F(v')dv'=F*F, (5.10) 

JO , J - oo 

an even function of v. Using this relation, (5.7) becomes 

W(v)=d 2 G(vW,b(v)+2F*F]. (5.11) 

Figures 7 and 8 show Q{v) and W(v) for a particular example. 

We wish to find the mean square fluctuations w 2 —w 2 in wit), the voltage after the modula- 
tion filter. It is obvious that 

and 

w=7jg(0) y (5.12) 

where g(v) is the amplitude response of G(v). For the case in hand, g(0)=0. 

The mean square voltage after the modulation filter is the integral of W(v) over all frequen- 
cies. The function F*F } which is related to the I F amplifier response I(v), is very slowly varying 
over the band of G(v). We assume that it may be replaced by its value at v m and removed 
from under the integral sign. Hence we have 

w*=a 2 G(0)4< 2 +2a 2 F*F\ Vm f G(v)dv. (5.13) 

We define the equivalent width of the function G(v) in the usual way, 

Ave= 2ck) f-J {v)dv - (5J4) 

Values of Av G for a number of filters are given by Colvin [1959]. The mean square fluctuations 
alter the modulation filter are then 

A^=^—w 2 ==4a 2 F*F\ Vm G(v m )Av G} (5.15) 

which is simply the area of the continuous spectrum illustrated b}^ figure 8. 
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The function F(y) is proportional to I(y) by (4.4) and (4.5), since Ii(v + p ) J r B(p— v ) is 
slowly varying where I(v) is significantly different from zero. 
The form of F(v) may then be rewritten as 



F(v)=A(T G /T +n)kT I(v)=ANkTJ(v), 



(5.16) 



where A is a constant of the instrument, T G is the equivalent temperature of the equivalent 
source at the receiver input terminals described in figure 3c and (3.32), and n accounts for con- 
stant noise. For the purpose at hand, it is sufficient to ascribe all this noise power density 
to an operating noise factor JV, in the usual way, as shown. We remark parenthetically the 
well-known fact that N depends on T G ; and the better the noise factor, the more sensitive is 
the dependence. From (5.16), it is evident that 

F*F\ Vm =(AkNT ) 2 I*I\ Vm . 
Writing Aw=(Aw 2 ) l/2 for simplicity, the output fluctuations are 

Aw=aAkNT [4l*Il m G(p m )Av G ) 1 ' 2 . 

Next we ask what output signal occurs when we intentionally introduce an equivalent 
temperature change AT G at the input to the receiver. We note that F, \p , and T[ will change 
by amounts AF=AkAT G I(v), A\p , and aA\p 0) respectively. If this change is carried out 
sinusoidally at the modulation frequency, assumed slow compared to the times over which the 
correlation functions and spectral densities are evaluated, a signal of amplitude w s will pass 
the modulation filter such that 



Ws=a\g(p m )\Axl/ =aAkAT G 



[G(v m )\ 1/2 f C _ 



I(v)dv 



This amplitude could be obtained more formally and elegantly perhaps, but surely less directly, 
by incorporating the effect of modulation into \p(r) and Q(v); but for this purpose the simpler 
argument will do. Following Colvin [1959], we may obtain the equivalent width of the func- 
tion /*/ as 



a ^ = 2tVo J.. pidv= 2M\o LL /(i H 



(5.17) 



using (5.16), (5.10), (5.8), (5.9), and (5.2). Values for a number of such receiver responses 
are given in figure 10, and are found to be closely related to the width of the IF response. 
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Figure 10. Bandwidths for several IF responses 
[after Colvin, 1959]. 



In terms of these quantities the ratio of mis noise fluctuations to signal amplitude is 

Aw/wM^NTjAT^iAv^A^Y^PIl, //*/| ) 1/2 . (5.18) 
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The third factor is very close to unity because v m is close to zero compared to the I V bandwidth, 
and we neglect its departure from unity. For signal to fluctuation ratio of unity, we find, 
dropping the subscript f/on AT G , 

AT=^2(Av (; /Ap I * i y /2 N7 () . (5.19) 

The analysis of an unmodulated radiometer employing d-c output and low-pass post-detection 
smoothing filter is contained in (5.13) and (5.12) and gives the same result as (5.19). The 
result essentially agrees with Dicke [1946]. 

So far the minimum detectable AT depends only upon two bandwidths and the system 
noise figure. It depends upon signal level only as this enters into the noise factor. In the 
spectrum of the fluctuations, the d-c terms have been suppressed by the use of modulation and 
the modulation filter, so the drift is eliminated to the extent that a-c amplifiers are more stable 
in gain than d-c amplifiers. Furthermore, gain variations at frequencies outside the bandwidth 
of the modulation filter do not affect the balance. In practical instruments, narrow bandwidth 
at v m may be obtained by twin-tee circuitry, tuned feedback amplification, or synchronous 
detection at this point. We note in passing that due to the spectral distribution shown in 
figure 7 having a triangular shape, there is in principal a reduction of background noise when 
v m -^B. In our case v m is relatively close to zero, and the effect is unimportant. In some 
applications sensitivity may be increased by raising the modulation frequency. 

Corresponding to a minimum detectable temperature change is a minimum detectable 
change in available noise power given by 

Al> mh =kATAv. (5.20) 

Equations (5.19) and (5.20) were used to compute the sensitivity of a receiver with a rectangu- 
lar pair reception filter and an ideal low-pass smoothing filter. The computed figures agree 
well with experimental values, indicating that the assumptions employed in this section are 
reasonable. 

It is not too difficult, following Lawson and Uhlenbeck [1950], to continue the analysis 
through the third detector and to examine the various possible combinations such as quadratic- 
quadratic, quadratic-linear, linear-linear, and linear-quadratic. Only qualitative results 
will be noted here, however. For the linear-linear and quadratic-linear combinations, the over- 
all rms fluctuation-to-signal dependence is of the form of (5.18) except with Av H) the width of 
the low-pass filter, replacing Av G , and thus enabling a narrower post-detection bandwidth. 
For the quadratic-quadratic and linear-quadratic cases, the overall rms fluctuation to signal 
dependence is as 

Aw^ /MW (Av G Av H )^ 
w' s ^KaTJ " Av^j 

The equivalent post-detection band is the geometric mean of the G and H filter widths. The 
criterion, rms fluctuation equal to signal, still gives an expression of the form of (5.19) for mini- 
mum detectable AT. 

6. Balancing Procedures 

Modulation of the elements p i3 - and q^ of the cascading matrices P and Q by any function 
of time produces a related time dependence of the power P in (4.7), through the two terms in 
which \a\ 2 and \c\ 2 appear. For the term in \a\ 2 the major time dependence is in the factor 
l/|pn| 2 . The quantity \T$i\ of (4.7) is slightly time dependent through $h (an element of the 
scattering matrix corresponding to the cascading matrix P) and through Y[ and Y' 2 which 
(Miter in T 3 i. However, the time dependence contributed by these factors may be verified to be 
at greatest of the second order in small quantities (that is, small scattering matrix elements or 
reflection coefficients). For the purposes of this work, this is assumed negligible compared to 
the zeroth-order time dependence of l/|^n| 2 . 
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In (4.7) we may then set the fundamental modulation frequency component of the power 
P to a minimum given by (5.19) and (5.20). The radiometer is used as in figure 11a with sources 
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Figure 11. Balancing procedure, Example A. 
Substitution of standard 7' 3 , Tz for unknown T\, T u 



T u I\ and T 2 , T 2 of (3.38) and (3.39) on arms 1 and 2 respectively. We take T l = T 2 =T Q 
for the temperature of networks P and Q. We designate the fundamental of the time-depend- 
ent factors by the subscript m. If the modulation responses of P and Q are carefully assured 
to be opposite in phase at the fundamental modulation frequency, balance requires that 

/ i \ IT^i-lr^iT.-T,) / i \ ]r 3 2l 2 (i-lr 2 i 2 )(r 2 -r ) 

\\pn\ 2 L li-r^l 2 VlaiilV* |i-r 2 ^| 2 {bA) 

We have recognized that the fundamental component of the reciprocal of a function of time 
is not in general equal to the reciprocal of the fundamental component of the function. Equa- 
tion (6.1) is important as the balance equation of the instrument, from which measurements 
and errors for any method of operation may be obtained. 

Example A. If we substitute a source T 3 , T 3 of absolutely known value, here supposed to 
be the standard, on arm 1 as in figure lib, and rebalance by adjusting only p n to a changed 
value p c n, we have the new balance equation 



\|pfi| 2 A 



| 2 (1- 



|r 3 | 2 )(r 3 - 



■To). 



i-r 3 ali : 



Alfliil 2 ! 



ir 32 l 2 (i-lr 2 ] 2 )(r 2 -r ) 
|i-iW 



(6.2) 



Note that IVh has gone into T 3 s\i because of the source interchange and the readjustment 
of the network P. The quantity |T 3 i| 2 depends on the reflection coefficient placed on arm 1, 
from (3.25) and (3.28) through (3.31). The dependence stems from terms only of second order 
and higher in small quantities (that is, reflection coefficients or small scattering matrix ele- 
ments). In practice this dependence may be reduced to terms of higher order if isolators, which 
tend to keep T[ constant although r x changes, are included in arm 1. Similar remarks apply 
to |T 32 | 2 . 

The source T 2} T 2 serves only as a reference comparison, and its characteristics as well as 
the characteristics of transmission in arm 2 are eliminated between (6.1) and (6.2). We obtain 
for source 1 in terms of source 3, 



(i-|r 1 | 2 )(r 1 -r )=^(x?/x 2 )(i-|r3[ 2 )(r 3 -r„), 

where we have abbreviated X 2 =(l/[p n | 2 ) m and 

^=!i-r 1 sj 1 | 8 /|i-r 3 sU| 8 . 



(6.3) 



(6.4) 



The error factor -n A requires knowledge of only the two additional quantities sh and s\i, both 
nearly equal and easily measured. If the matching is good, r) A is very nearly unity. This is 
the mode of operation presently in use at the National Bureau of Standards. 
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As a practical point, it is worth Doting that the choice of the transmission characteristics 
of the three-port junction may be chosen to advantage in special cases. For use as in Example 
A, we may choose 4 a directional coupler rather than a hybrid tee and secure the required throt- 
tling of T 2 hy the coupling factor instead of by an attenuator in network Q. This design avoids 
an approximate :; db power loss of the low-temperature source in the junction and improves 
the sensitivity as described below. The present radiometer analysis is general enough to 
apply to this design also. 

When only comparison of sources at high temperature (~ 11,000 °K) is required, sensitivity 
of comparison is increased by the ability to keep the attenuation ratio in network P about unity 
instead of throttling the high-temperature sources down to the level of the low-temperature 
(-1,000 °K) hot-body standard source. Then a change AT X or AT 2 in the sources is trans- 
mitted to the receiver unattenuated, and this quantity inust be equal to or greater than the 
minimum detectable change given by (5.19). Otherwise only the attenuated change reaches 
the receiver, and correspondingly greater change at the sources is required to produce the 
minimum detectable change at the receiver. 

Example B. The constants of the instrument, T 31 , T 32 , and q u , may be eliminated by 
interchanging sources 1 and 2 in (6.1) and adjusting p n to a changed value p c n for rebalance, 
as in figure 12, which also changes s\ t to s\i. Elimination of the above constants between 
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(6.1) and the corresponding equation describing the interchange gives for source 1 in terms 
of source 2, 

(i-|r 1 | 2 )(T 1 -T )=^(x c /x)(i-|r 2 |2)(r 2 -r ), (6.5) 

where 

^Hl-IVhl ll-IY^I/jl-IVlfl |l-r 2 ^|. (6.6) 

Again, y\ B is very close to unity for good matching. In this method, evaluation of n B requires 
knowledge of the additional quantity s 2 n , also easily measured. This method is best used for 
sources of similar T, T. For, if the temperatures are too dissimilar, it requires excessive throt- 
tling in both arms, resulting in poor sensitivity. But the more similar the source reflection 
coefficients, the closer rj B becomes to unity. 

Example C. An alternative method of operation may be convenient when many unknown 
sources are to be compared in succession with the standard, T 3 , T 3 . As a preliminary adjust- 
ment, trial sources 1 and 2 of similar T and r are repeatedly interchanged on arms 1 and 2, 
and p n , <Zn> an d T 2 are adjusted until further interchange does not disturb the balance. Then 
in (6.1) A=A C and sli=sli, and we may solve for the instrumental constant 



: \kiil 2 A 



li-r^lli-ivi; 



(6.7) 
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The standard source T Zj T 3 may then be permanently installed on arm 2 as in figure 13. Then 
for any source T, T installed on arm 1 and balanced to the standard without further adjust- 
ment in arm 2, but only by adjustment of p c n to a new value p n (and s\i to Sn and X c to X), 
we have from the form of (6.5) corresponding to this situation and from (6.7), 



(i-|r| 2 )(r-r )=, c (x?/x 2 )(i-|r 3 | 2 )(r 3 -r ), 



where 



Vc= 



li-ri&Hi-iv&l |i-r«},| 



i— rx*if l— rasij l— ras?, 



(6.8) 



(6.9) 



Again, n c is very close to unity for well matched units. 



t, r 




P 




S 




Q 




V 3 


~ll,000°K 




4TTEN-I/ 


1 








ATTEN ~l 




-1,000° K 



Figure 13. Balancing procedure, Example C. 

Comparison of unknown T, r with T3, T3, having previously cal- 

brated the instrumental constants. 
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These examples serve to illustrate the numerous methods of possible operation, not all of 
which are independent. For example, it is stated without proof that Example C reduces to 
Example B for T, Y=T U T u and T 3 , T 3 =T 2 , T 2 . Also Example C reduces to Example A if 
T 3 , T 3 replaces T, T in the left-hand member of (6.8) and T 2 , T 2 replaces T 3 , T 3 in the right- 
hand member. 

7. Typical Error Analysis 

So far we have not assumed any specific configuration for the arms of the radiometer. 
In order to demonstrate that a ratio such as (l/\p c n \ 2 )mK^/\pn\ 2 )m = K/^ 2 is an attenuation 
ratio, and to evaluate error factors numerically, we choose a particular case where one arm 
consists of an attenuator and a modulator. Addition of isolators before or after this combina- 
tion will not affect the result because of the properties of isolators. 

We write a cascading matrix of the type denned by (3.18) for an attenuator with atten- 
uation IZI 2 , 



l -r z2 

Vn I T n T i2 _ 



(7.1) 



where we have used the reciprocal nature of the attenuator. 

We wish to write the cascading matrix M for the general case of a nonreciprocal modulator 
in such a way that we may easily find the fundamental component of l/|^n| 2 . The scattering 
matrix is 



S m = 



r«i(0 

_m 1 (t) 



rn 2 (t) 

r„»(0J 



(7.2) 



We assume the elements are constant functions of frequency over the RF pass band, and we 
display the slow time dependence at the modulation frequency. This time dependence in 
practice is, unfortunately, generally different in form and magnitude for each element. From 
(7.1), (7.2), and (3.19) we may form the total cascading matrix, find the element p n = (LM) n , 
and finally find 

1 _ ZraiOO 



Pu 
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If we take the two time dependences to be real and separable, of the form ?n i (t) = m 1 f(t) and 
T ml (t) = T ml g(t), we see that 

' NZIVil 2 (L/(*)] 2 +O(r 2 )<7(0[/(*)] 2 }, 



bn 12 



where the terms through second order in small quantities (the reflection coefficients) times 
g(t) are displayed. When the fundamental component X 2 of each term is taken, to a very high 
degree of approximation, 

K/V=\lc\ 2 /\l\ 2 

because of almost complete cancellation of the second-order terms. This cancellation is the 
more complete the more constant T I2 is with attenuator setting. Introducing a and 7, the 
attenuator readings in decibels corresponding to \l\ 2 and \l c \ 2 , respectively, 



IOTI 2 =io- 



(7.3) 



Numerical evaluation of our results will permit examination of the relative importance 
of each source of error. First, the attenuator calibration may be certified to ±0.01 db. Sec- 
ond, the error in measurement of the standard temperature T 3 may be ±5 °K or ±0.022 db. 
The output fluctuations given by (5.19) can be made negligible compared to other errors. 
Third, assuming the numerical magnitudes |r| = |r 1 | = |r 2 | = |r 3 |=0.030, I^H^nH |s?_|=0.050 
(any differences may be accounted for more precisely in practice), and allowing for the most un- 
fortunate phases, we find from (6.4), (6.6), and (6.9) that the bounds on the 77's are \rj A \< 0.026 
db, \r] B \ < 0.026 db, and \rj c \< 0.052 db. If one did not know or did not wish to carry the 
reflection coefficients associated with each source, one could write (6.3), for example, as 



T.-T^n^WW.-T,), 



(7.4) 



where we take ^=^4(1 — |r 3 | 2 )/(l — jr^ 2 ). If this error limit were tolerable, (7.4) could be 
used for the radiometer. Imperfect matching may in practice contribute significant error, so 
that it is evidently desirable to match the system as well as possible. 

Table 1 gives a comparison of an early model (1958) NBS radiometer with a still earlier 
(1946) MIT radiometer and with a maser-TWT radiometer at Columbia. The measured out- 
put fluctuation of the NBS radiometer was ±0.2 °Iv, which agrees roughly with (5.19). This 
indicates that our assumptions are in some measure reasonable. 



Table 1. Calculated performance of three radiometers 



Performance characteristic 


MIT 


NBS 


Columbia « 


Receiver noise figure. 


25 
8 

0.4 
0.95 

1.0x10-9 


3.2 

10 
0.2 
0.08 
0. 96X10-10 




IF bandwidth Ai>i*i 

Low-pass filter taw 2irAuo 

Calculated AT 

Calculated minimum detectable power change.. 


Mc/s__ 

sec -1 .. 

°K__ 

__ergs/sec__ 


5.5 

0.02 
0.15X10-1" 



' Drake etal. [1958]. 



8. Summary 

It is possible to derive a very general expression for the voltage wave emergent from one 
port of a multiport microwave junction, which itself contains sources, and to the arms of which 
are connected general microwave two-ports themselves terminated by active sources. This 
expression may be specialized to a class of microwave radiometers consisting of a three-port 
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junction, to the arms of which are connected, through general microwave two-ports used for 
measurement, respectively a microwave source to be measured, a standard source in terms of 
which the measurement is made, and a receiver. The analysis is sufficiently complete to 
identify clearly many conditions which bear on the accuracy of measurement of the unknown 
source, and to quantitatively evaluate them, or to keep them below some tolerable contribution. 

In particular, noise arising from lossy elements anywhere in the microwave network can 
be explicitly taken into account. The familiar result obtains, namely that the noise associated 
with the excess of the source temperature over the radiometer temperature is observed ; but 
this result is now subject to the important condition that all lossy elements of the radiometer 
proper be at the same temperature. Adequate isolation among arms is required to avoid 
contributions from arms other than the desired one. Departures from ideal impedance match 
in the microwave arms require corrections, which are here evaluated to second order in small 
quantities ; this analysis constitutes the main quantitative result of the paper. For radiometers 
designed for rapid comparison of two sources, exactly opposite phase in the switching or modula- 
tion is necessary. Isolation against impedance changes in each arm becomes necessary for 
valid measurement. Correct treatment of the image response of the receiver is also necessary. 

An analysis of fluctuations extends previous results for radiometers which are designed 
for the continuous observation of a source, and which therefore have low pass filters after the 
second detector. The extension applies to radiometers which are designed for rapid com- 
parison of two sources by switching or modulation, and which therefore have band pass filters 
after the second detector. The sensitivity to source temperature change is dependent on the 
receiver noise factor, the limiting IF bandwidth, and the low frequency filter bandwidth in the 
same way as previously. The sensitivity of radiometers with a third detector operating on the 
signal at the switching or modulation frequency depends on the particular combination of 
linear or quadratic second and third detectors used. The spectral distribution of the fluctua- 
tions at various circuit points is useful as an aid to understanding and design. 

A general balance equation for the instrument, subject to explicit simplifying conditions, 
may be derived and specialized to three convenient modes of operation in order to account 
for constants of the instrument appropriate to the two input paths. These modes are (a) 
substitution of a standard source for an unknown source on a particular arm, (b) interchange 
of a standard source and an unknown source on the two arms, and (c) calibration of the instru- 
ment so as to allow immediate comparison of an unknown source with a permanently affixed 
standard. Measurement is accomplished in terms of only the s 2 i (or p n ) element of the equiv- 
alent scattering (or cascading) matrix of the measurement arm. 

Corrections to the measurement explicitly appear in terms of a relatively few pertinent 
reflection coefficients and scattering matrix elements. These expressions may be used either 
to correct the observations or to keep the errors within acceptable limits. Typically, reflection 
coefficient magnitudes of a few hundredths require correction of a few hundredths db. An 
important result is that mismatch of the standard and unknown sources to the waveguide 
transmission line requires correction only as the difference of the absolute squares of their 
reflection coefficients, and so may often reduce the severity of error if the sources are similarly 
mismatched. 

Radiometers in use at the National Bureau of Standards, to which this analysis is applic- 
able, have been previously described by Estin, Trembath, Wells, and Daywitt [1960], and by 
Wells, Daywitt, and Miller [1962]. 



Early work on NBS microwave noise standards was carried out by J. J. Freeman and 
C. R. Greenhow. Experimental work was carried out by M. M. Anderson. Acknowledge- 
ment is due W. R. Atkinson, R. W. Beatty, E. L. Crow, A. J. Estin, and J. S. Wells for their 
suggestions. 
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